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Abstract The catastrophic senescence of the Pacific salmon is among the initial tests used 
to validate the Penna aging model. Based on the mutation accumulation theory, the sudden 
decrease in fitness following reproduction may be solely attributed to the semelparity of the 
species. In this work, we report other consequences of mutation accumulation. Contrary to 
earlier findings, such dramatic manifestation of aging depends not only on the choice of 
breeding strategy but also on the value of the reproduction age, R, and the mutation thresh- 
old, T. Senescence is catastrophic when T < R. As the organism's tolerance for harmful 
genetic mutations increases, the aging process becomes more gradual. We observe senes- 
cence that is threshold dependent whenever T > R. That is, the sudden drop in survival rate 
occurs at age equal to the mutation threshold value. 
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1 Introduction 

Semelparous is the term given to species that undergo genetically programmed degeneration 
following procreation [ 1 1. They breed once at the same age, usually producing plenty of off- 
spring, and die shortly after. These semelparous organisms include annual plants and many 
species of animals, particularly in the class Insecta |2]. It is generally believed that these 
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organisms allocate a huge amount of energy towards reproduction, thus making death in- 
evitable 1 3 1. Unlike in populations that breed repeatedly (iteroparous), senescence associated 
with semelparity or one-time reproduction is abrupt, and is often described as catastrophic. 

The life of a Pacific salmon, for example, traces a path from freshwater where they are 
born, to the ocean where they mature, and then back to freshwater where they spawn (4). 
Compared to other species which repeat this cycle, the Pacific salmon is only capable of 
one full roundtrip. They die in freshwater just after producing offspring. One probable cause 
for such behavior is starvation. Literature indicates that the salmon sometimes travels back 
more than 1200 kilometers to freshwater, generally without eating [5 |. However, according 
to the mutation accumulation theory, the catastrophic senescence of the Pacific salmon is 
due solely to the semelparity of the species |6 |. That is, the drastic lowering of the salmon's 
fitness (to the point of death) immediately following procreation is completely explained by 
reproduction occurring only once at a fixed age. 

Another study |7 1 demonstrated such dramatic aging manifestation utilizing a model that 
does not account for the effects of harmful mutations. Instead, some form of antagonistic 
pleiotropy was employed, promoting an increase in births by either shortening individual 
lifespan or delaying reproductive maturity. 

This paper presents a further investigation on mutation accumulation and semelparity 
within the framework of the Penna model [8|, a simple, easily extensible tool for predicting 
many of the aging-related features found in biological systems [9|. Our results show aging 
trends other than the previously observed catastrophe. 



2 Methodology 

We implement the Penna model described in 15 1 101 . Individual characteristics are stored 
in 32-bit long genomes. For every year in the individual's life, one bit in the genome is 
read. Thus, only those genes located at bits less than or equal to the individual's current 
age are considered active. Zeroes correspond to healthy genes while ones indicate harmful 
mutations or diseases. The mutation threshold, T, defines the upper limit on the number of 
mutations individuals can undergo. An individual suffers a genetic death when the number 
of active mutations reaches this value. Reproductive maturity is achieved at age R. The num- 
ber of offspring produced by a parent at each time of reproduction is given by the birth rate, 
B. Genetic traits, both active and inactive, are passed on from parent to offspring. During re- 
production, M new harmful, randomly located irreversible mutations are introduced into the 
newborn's genome. Usually, the mutation rate, M, is set to one [?]. Environmental restric- 
tions are taken into account by imposing an age-independent Verhulst factor, V, = \ —N t /K. 
This corresponds to the probability that an individual survives given the current population 
size, TV,, and the carrying capacity, K. In simulations, the carrying capacity is commonly set 
to ten times the value of the initial population. 

Simulation details A population of 20000 is propagated in an environment that can support 
a maximum of 200000 individuals. Initially, all genomes are clean (no defective genes) 
and all individuals are at age zero. To consider semelparity, reproduction age is fixed at 
R. Because such species are generally known to produce plenty of offspring, the original 
paper utilized high birth rates. However, since the catastrophic behavior was shown to be 
independent of the birth rate | 6 |, we limit its value to ten {B = 10). Furthermore, we extend 
the demonstration to higher thresholds. 



3 



(a) (b) 




1 1 1 1 

8 16 24 32 

Age 

Fig. 1 Survival rates as a function of age in the case of reproduction happening only at R = 17. The sharp 
decline in organism fitness occurs at R = 18. 

The simulation runs for 3000 iterations. Simulated populations have already achieved 
steady state before this time. Demographic statistics such as the fraction of defective genes, 
Verhulst-related mortality and survival rates as functions of age are obtained for steady state 
populations associated with different parameter values. To minimize fluctuations, we con- 
sider mean population sizes of the last 300 iterations. Data hereafter are representative of 
ten independent simulation runs. 



3 Results and discussion 

Aging is said to be catastrophic for populations that maintain a high level of fitness and 
condition until shortly before death 1111 . Losses in functional capabilities directly affect 
an organism's ability to withstand a given environment. Thus, here, we relate fitness to the 
survival rate - the ratio between the population with age a at time t and the population with 
age a — 1 of the previous iteration, N a (t) /N a -\{t — 1). 



3.1 Catastrophic senescence 

The survival rate of semelparous species usually goes down, from some relatively high value 
to zero, at the age immediately following reproduction. For populations associated with 
R = 17 and T < R, the sudden decrease in organism fitness happens at R = 18 regardless of 
the value of T (Fig. [TJ. 
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Fig. 2 Genetic profile of steady state populations associated with R = 
higher bits, starting at 18. 



17. Harmful mutations accumulate at 



Catastrophic senescence is closely associated with the accumulation of bad mutations 
at the end part of the genome |12|. The life cycle of semelparous species can be viewed 
as a two-phase process driven by different mechanisms |13j . It starts with a juvenile non- 
reproductive phase (age < R), followed by the adult phase which is characterized generally 
by reproductive maturity. Genes switched on during the juvenile phase are called house- 
keeping genes, while those activated in the adult phase are referred to as death genes [ 14 1. 
The concept of mutation accumulation introduces strong selection pressure during the non- 
reproductive phase in order to regulate harmful mutations that are passed on from parent to 
offspring [5|. Genetic housekeeping works to minimize harmful mutations expressed prior 
to reproduction. Selection, thus, acts to keep the first part of the genome clean. When T = 1, 
only individuals with no expressed mutations at bits less than or equal to R are allowed to 
procreate. The model is designed, however, in such a way that inactive genes are also inher- 
ited. The selection process brought forth by the dynamics of mutation accumulation does 
not screen bad mutations that are activated only after reproduction. As a result, diseases pile 
up at the higher bits. Fig.[2]shows the corresponding genetic profile of the steady state popu- 
lations in Fig.[T] The fraction of deleterious mutations occurring at bits greater than R is 1.0. 
Moreover, we see that when T < R the survival rates fall off at the age when the fraction of 
defective genes is equal to 1.0 (cf. Figs. [T] and Thus, the catastrophic loss in organism 
fitness is a consequence of the high density of bad mutations switched on immediately after 
reproduction. 

The strength of selection diminishes with increasing mutation threshold. At higher T, 
more mutations are preserved in the gene pool. Such genes are believed to be essential for the 
organism's growth. It was previously demonstrated [6] that although higher threshold values 
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Fig. 3 Survival rates (solid lines) and genetic profile (or the fraction of defective genes - circles) of semel- 
parous populations associated with R = 27. Senescence is more gradual at higher threshold - (a) T = 21, (b) 
T = 23, (c) T = 25, and (d) T = 27. 



yield larger populations, a different T does not alter the survival rate. On the contrary, Fig. [3] 
shows that as T is increased further, the manifestation of aging becomes more gradual. The 
abrupt decline in organism fitness still occurs at R + 1. However, a fraction of the population 
is observed to survive after procreation. 



3.2 Threshold dependence 

When T > R, selection is very weak. Genetic housekeeping is perfomed solely by the Ver- 
hulst factor because threshold deaths do not occur until age T. In contrast to mutation accu- 
mulation, the Verhulst factor acts at random. It is incapable of discriminating between good 
and bad genomes, killing individuals regardless of fitness. It does very little in preventing 
the spread of diseases among the bits. Unlike in the previous case where bad mutations are 
more concentrated at higher bits, we observe a high density of deleterious genes occurring 
even at lower bits. In Fig. [4] the fraction of defective genes is approximately 1.0 for all bits. 
For these populations, the drop in the survival rate is manifested as soon as organisms reach 
age = T. Senescence, in this case, is threshold dependent. 

Plots in Fig.|4]also register a local minimum for the survival rates at R + 1 . The frequency 
of such dips increases with lower R and higher T. Fig. [5] presents the survival rates of pop- 
ulations associated with R = 1,T = 16,24 and 32. The corresponding bottom plots show 
Verhulst-related death rates. Notice that the dips happen at ages when the Verhulst mortality 
rate is highest. Recall that in our simulations, we begin with a set of individuals character- 
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Fig. 4 Behavior of the survival rates when T > R. Here, R = 17: (a) T = 20, (b) T = 24, and (c) T = 28. 
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Fig. 5 Survival rate (top) and Verhulst mortality rate (bottom) of steady state populations associated with 
R = 7: (a) T = 16, (b) T = 24 and (c) T = 32. 
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ized by the same age and genome. The resulting population is, thus, structured in such a 
way that those with ages separated by a factor R are clustered together. The Penna model 
follows Gompertz law 1 9 , 15 ] which implies that most individuals are newborns. Because the 
Verhulst killing effect increases with total population size, it is felt most by those grouped 
with age 1. The observed periodicity, hence, arises from our choice of initialization. This is 
not seen in other Penna model implementations whose initial population is characterized by 
random ages and genomes. 

We lose the dips by limiting Verhulst deaths to newborns only - VB implementation 
[ 17 1, as in Fig. [6] Note, however, that in our demonstration, we had to use a lower birth rate, 
B = 3, in order to avoid extinction via overpopulation. Higher B values introduce sudden 
fluctuations in the total population (beyond the set carrying capacity) which kills all new- 
borns and, eventually, the entire population. The maximum lifespan is generally longer for 
the VB population |16|. For the above presented semelparous cases, this is evident only at 
larger R values (cf. Figs.[4][5l and [6j. 

4 Conclusion 

The theory of mutation accumulation depicts senescence as a consequence of increased 
mutation load occurring at later ages [ 18 1. Within the framework of the Penna model, this is 
seen in the high concentration of bad genes located at bits greater than R. An increase in the 
tolerance for harmful mutations weakens selection and increases the probability of finding 
a deleterious mutation at a lower bit. When T > R, this probability is equivalent to unity, 
for all bits. Ignoring Verhulst effects, the corresponding survival rates are henceforth flat at 
1.0 from ages 1 to T — 1. Such lack of age-specific variation in fitness and mortality rate has 
been reportedly observed in coldwater fish, bivalves, turtles, whales, and naked mole-rats 
fT9ll20T 

Aging manifests differently depending on the choice of reproductive strategy. Senes- 
cence associated with semelparity is generally catastrophic, while that with iteroparity (or 
multiple breeding) is gradual. This characteristic makes organisms that exhibit interspecific 
variations in parity very useful in evolutionary studies [(1). The catastrophic behavior of the 
Pacific salmon has been modeled extensively with tools that are based on different aging 
theories. Previous works [ 6 , 7 ] suggest that the one-time reproductive strategy is the only in- 
gredient needed to explain the substantial loss in functional abilities immediately following 
procreation. Utilizing a model that is based on the mutation accumulation theory, we were 
able to find other limiting factors that yield catastrophic senescence. When T < R, the high 
density of bad mutations occurring at higher bits causes an abrupt decline in the survival 
rate at age R + 1 . On the other hand, when T > R, aging is catastrophic at age T (threshold 
dependent). The major drop in fitness happens at age T, when threshold deaths begin to take 
place. 

References 

1 . Crespi BJ and Teo R, Comparative phylogenetic analysis of the evolution of semelparity and life history 
in salmonid fishes, Evolution, 56, 1008-1020 (2002) 

2. Coleman BD and Ying-Hen Hsieh, Theory of the dependence of population levels on environmental his- 
tory of semelparous species with short reproductive seasons, Proc. Natl. Acas. Sci. USA, 76, 5407-5410 
(1979) 

3. Smith FA and Charnov EL, Fitness trade-offs select semelparous reproduction in an extreme environment, 
Evolutionary Ecology Research, 3, 595-602 (2001) 



8 




(b) 



0.5 




16 
Age 



24 



32 



T=20 
T=24 
T=28 
T=32 



Fig. 6 Survival rate resulting from a VB implementation: (a) R = 7 and (b) R - 
rate was reduced significantly to avoid extinction via overpopulation (B = 3). 
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